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^ ' Abstract 

(<-^ ' The valuation of counterparty risk for single name credit derivatives 

pf-\ , requires the computa- tion of joint distributions of default times of two 

f^ ' default-prone entities. For a Merton-type model, we derive some formulas 

for these joint distribu- tions. As an application, closed formulas for coun- 
terparty risk on a CDS or for a first-to-default swap on two underlyings 
are obtained. 



1 Introduction 



Default Risk has become one of the key issues of contemporary finance. In 
C^ I practice, defauh risk computations are directly involved in three, closely con- 

nected, but still separated areas: first, the pricing of single and multiname 
credit derivatives (ABS, CDS, CDOs...), second, risk valuation for regulatory, 
risk management and economic capital valuation purposes, third, rating as- 
signements. The present article is focused mainly on the first two aspects of 
risk valuation. Notice however that the next generation of rating methodologies 
will most probably have to include advanced correlation tools such as the ones 
considered here. 

Although deriving new formulas for credit derivative transactions on two 
assets, our main interest will be on counterparty risk. Once again, this is a 
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central argument, both for regulatory purposes (counterparty risk was already 
a key issue in the Basel 1 agreements, but its relevance and complexity was 
emphasized by the later improvements of the set of rules, see e.g. [I1IH1I2]) and 
also for practical risk management. 

Concretely, the 2007 "subprime" credit crisis has emphasized how a mis- 
pricing of counterparty risk could be dangerous. Before the crisis, financial 
institutions engaged heavily in credit derivatives transactions for hedging pur- 
poses. One of the main side effects of the crisis has been to reveal that their 
counterparties were not reliable: monoline insurers (who act traditionaly as 
guarantees for municipalities-type bond emissions but engaged more recently in 
the ABS business) faced downgrades and/or bankruptcy, whereas some compa- 
nies offering CDS protection have also been on the wedge of bankruptcy due 
to spread widenings and collateral agreements. In spite of these facts, in the 
present state of the art quantitative counterparty risk assessment is still largely 
in infancy. 

We address here one of the fundamental questions in the field, namely we 
introduce a Mcrton-typc model and derive formulas for the joint distributions 
of the default events -and other relevant joint probabilities- of two default-prone 
entities. Closed formulas for counterparty risk on a CDS transaction or for a 
first-to-default swap on two underlyings follow. 

Notice that we refrain to seek for the outmost generality or for a complete list 
of pricing formulas that could be obtained using the methods developed in the 
present article. Most credit derivatives on two underlyings with time-dependent 
payoffs can actually be priced using our results. The exercise of deriving the 
corresponding formulas is left to the interested reader. 

2 Modeling two dimensional default risk. 

Recall that we are interested in modeling the counterparty risk on a credit risk 
transaction. This involves (in general and depending of the particular features of 
the transaction) the computation of various joint laws. For example, for a vanilla 
credit default swap, the computation relies on the distribution of the default 
time of the counterparty and the conditional distribution (w.r. to the default 
of the counterparty) of the values of the CDS contract. The same methodology 
can be applied to the valuation of first-to-default contracts on two underlyings 
when the payoff is time-dependent (see section 4 of the present article; the 
time-independent case has been addressed in p[T71fT3]). 

We follow Merton's structural approach to risk, where the default of a com- 
pany is triggered by its firm value falling below a threshold (see [H] or e.g. [H 
Chap. 3]). It is well-known that the original Merton methodology, where the 
threshold is determined by the firm's liabilities, as well as its natural refinements 
such as the Black-Cox first passage time formulas j3], give qualitatively good 
but numerically poor results. For example, the structural method will rank cor- 
rectly the risks of two companies but will not be able to predict the spreads of 
the corresponding CDS. For pricing purposes, the correct way to use the Merton 
and Black-Cox models is to calibrate their parameters on the market prices of 
the securities issued by the companies. Doing so insures that the models price 
correctly the risk. This is the ground for their various uses, cither in credit risk 
measurement (think to the Basel 2 Vasicck-type large pool formulas [HI [5]), 



either for pricing purposes (think to the one-factor Gaussian copula model for 
CDO tranches [H]). 

Notice that the computations we are interested in require to go beyond 
the existing formulas for two-dimensional credit risk in the structural approach 
[8l[l7l[T3j. Indeed, the formulas in these articles rely mainly on the computation 
of the probability that one (and only one) out of two firms defaults or that the 
two firms default before a given maturity, whereas pricing counterparty risk for 
a CDS transaction requires (among others) the exact knowledge of default times 
distributions. 

We have tried to keep some balance between generality and the choice of too 
restrictive assumptions. Notice that the model may be simplified by selecting a 
suitable subset of meaningful parameters, depending on the particular situations 
and/or the available set of data. In the end, the model we use is as follows. 

We label 1 and 2 the two default-prone entities. In view of applications to 
counterparty risk, the counterparty will be labelled 2. We write ri (resp. T2) for 
the random time when the first (resp. second) entity defaults. Since we work 
in the structural model, defaults arc triggered when two log-normal processes 
Vi and V2 associated respectively to the first and second entity fall below a 
barrier process. As in Black-Cox 0, we assume that the barrier is described by 
a time-dependent deterministic process Vi{t) := Kie'^^^. At last, we assume that 
the risk-neutral dynamics of the processes Vi and V2 are given by: 

^^ = (r-h)dt + a,dB,it), 

where r is the constant short-term interest rate and i?i, i = 1, 2 arc two Brown- 
ian motions. The coefficient ki is a payout ratio representing net payouts/inflows 
by the firm, sec e.g. the account of Merton's approach in [H Chap. 2]. Notice 
that we use a risk-neutral dynamics since we are interested in pricing formulas, 
but the same construction holds with a real-life dynamics (just replace r — fc^ 
by the drift under the historical probability). 

At last, we assume that the Brownian motions Bi and B2 are correlated: 
Cov{Bi(t), B2{t)) ~ gt. This is an important assumption that reflects the fact 
that the counterparty of a CDS transaction is usually a well-rated company 
that will not default excepted in a very bad macroeconomic environment. In 
particular, if the counterparty defaults, one may expect CDS spread on other 
entities to widen considerably simultaneously. This phenomenon is accounted 
for by a positive correlation coefhcient. It is the main reason, together with 
the dynamical aspects of the problem, why simplistic counterparty risk models 
fail, both theoretically and empirically, to account for counterparty-driven loss 
expectations on derivative contracts. 

3 Counterparty risk 

In this section, we derive a closed formula for the counterparty default leg of a 
CDS contract, that is the present value (PV) of the expected losses on a vanilla 
single-name credit derivative transaction due to the default of the protection 
seller. Notice that the notion of counterparty default leg should not be confused 
with the usual notion of default leg for a CDS contract, that is, the present value 
of future payments by the protection seller. 



We follow standard market practices. Namely, we assume that the CDS 
contract has a notional value C and that, if the underlying entity defaults, the 
buyer of protection receives (1 — i?„) • C from the seller of protection, where Ru 
stands for the recovery rate of the underlying entity. Similarly, if the counter- 
party (the protection seller) defaults at t on the CDS contract, we assume that 
the buyer of protection receives (1 — R^ ■ V^ ^ where V^^ = sup(0, Vt) stands for 
the positive part of the market value of the CDS contract at t, and where Re 
stands for the counterparty recovery rate. The present moment is normalized 
to i = in what follows, it may or may not coincide with the inception of the 
CDS contract. 

3.1 Counterparty default leg 

Using the joint law of the default time of the counterparty and of the underlying 
entity for the CDS, we compute the default probability of the underlying entity 
conditional to the default of the counterparty, and obtain a closed formula for 
the counterparty default leg. 

The coupon payments rate (CDS spread) by the protection buyer is written s. 
We assume continuous payments of the fees: in practice, fees are payed quarterly 
(and sometimes semi-annually), but due to the payments- in-arrears conventions, 
the continuous payments assumption is a reasonable one. It is usually defined, 
at inception of the contract, by equalizing the fee and default legs of the CDS 
contract (see Section 13.31 Some empirical adjustment is sometimes done for 
the counterparty default risk, taking into account advanced risk management 
parameters such as netting agreements. Our computation of the counterparty 
default leg leads to a new, sounder, methodology, to assess the fair value of this 
adjustment. 

First, the theoretical counterparty default leg Dc is given by 

D, = (1 - R,) ■ E[e--^- ■ sup(0,p(V^i(T2),r2))l,,<(TAri)] 

where p(Vi(t2),T2) is the market price of the CDS contract at t = r2 (when 
Ti > T2). The value p{Vi{t2)tT2) is obtained as the difference between the 
default and fee legs, that is: 

p(^i(r2),T2) = A(V^i(r2),T2) - — • £;[(1 - e-'-«^^-^)--^))l.,>.,|^.,] 

r ~ 

where we write Tt, as usual, for the natural filtration of the probability space 
underlying the two Brownian motions Bi(t),B2{t). 

Notice that the term T A ti can be, in most situations, safely replaced by T 
in the previous expansion. This is because, for standard values for the spreads 
and implied default probabilities in single- name default risk computations, the 
computation of the fee leg conditional to the hypothesis that no default occurs 
is a good first-order approximation to the unconditional fee leg. The following 
computations could be simplified accordingly. 

Now, the value of -D;(V^i(t2), T2) is given by: 

A(Vi(r2),r2) = E[C{l-Ru)e-'^^'-^'-hr,<r,<T\Tr,] 



Theorem 3.1 The counterparty default leg Dc is given by: 

D, = C(l - Re) 
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The Theorem follows by standard computations of integrals over Gaussian 
densities that are sketched in the Appendix 1. 

3.2 Explicit formulas 

Recah that n = mi{t,V,(t) < X.eT'*}. The condition V^(t) < Kie'^-* can be 
rewritten: W,{t) > yj, where W,{t) = H^^^) - Hi^) and yl = In ^,(0) - 
In Ki. Equivalently, Wi{t) is the diffusion process: 

dW,{t) = -iy,t~a^dB,{t), 

with Wi(0) = and I'i := r — ki — 7,; — \(j'j. 
Let us define Z(t) by: 



Z(t) = (Zi(i),Z2 (<))* = 



We get: 
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dZi{t) ^(j)idt + dXi{t), dZ2{t) ^ (j)2dt + dX2{t), 
where X(t) is a standard planar Brownian motion and 

via2 - V2crig , V2 

In particular, Z(t) is a 2-dim. Brownian motion with drift and the barrier 
conditions Viit) ~ Vi{t) now read: Z2{t) = and -^/l — g'^Zi{t) + gZ2{t) ~ 0. 



We want to compute the default leg Dc, that is, equivalently: 

E[h{Bi{T2),T2)lr2<{TA.Ti)] = -E[/l(Zi (T2) , T2) 1^2<(TAri)] 

rT 1.+00 
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where h{x,t) := 

y r \ \/T ~ t \/T -t J 



^^^ ■= a^\i^it + crix + lnVi{0)-\nKi) &ndh{z,t) ^ h C^'"'* + ^1 - p^ ■ z,t). 
Applying the Girsanov theorem, (Zi(t), Z2{t))* is a classical Brownian motion 
for the probability law Q: 

Let roe^^o := Zi(0) + 2^2(0) = '^»'"^~f»'"^ + i ^. 
Lemma 3.2 We have, for (a, 0) € R^ s.t a > : 

P(t2 e di,T2 = r2 Ari,Zi(T2) € da) = e'^i(''-''«™"(^''»-'*^''°''''^(^°)-^^^e-("'+''o)/2' 

E°° . mrOo (O.'To 

nsm I„^/a{—-)dadt. 
at 
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where ||(^|p := 0^ + 02? ^ •" arcsin(gi) + ^ and In-KJa ** i/ie modified Bessel 
function of index nil /a. 

The Lemma follows from Thm. in Appendix 2 (up to a straightforward 
adaptation since we consider here a polyhedral domain D := {(.t, y) S R^|y > 
0, yl — p^a; + gy > 0} with a horizontal lower -instead of an upper- boundary, 
so that the signs have to be changed accordingly in the formulas). 

Indeed, according to [S] (see also [131 pp.697]), we have, for (a, b) in a neigh- 
borhood of (a, 0) with & > 0: 

f{a,b,t)dadb^ Cl{Z{t) e {da,db),Ti At2 > t) 
at ^-^ a at 

T1=0 



where 9 := arctg{-), fi := \/a^+l^. Therefore, we have: 
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and the proof follows. 



Theorem 3.3 The counterparty default leg Dc of the CDS is given by: 



Dc = C{l-Rc) / / /i(^^f)e'^i(^-'-o™«(eo))-<^2r-osi„(eo)-J 
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Notice that although involving a double integral on Gaussian densities and 
Bessel function (which asymptotic behavior is well-understood, making numer- 
ical approximations easy and efficient), this kind of higher order integrals is 
highly familiar in physics (heat conduction in solids, cross-sections computa- 
tions, acoustics...) and can be handled efficiently. However these techniques are 
not, for the time being, of our competence, and we postpone to further work the 
numerical analysis of the problem. For the use of Bessel functions in physics, 
we refer to the classical treatises such as O El [16] . 

3.3 Fair price of a CDS 

In this section, we take advantage of our computation of the counterparty default 
leg to compute the fair price of a CDS contract at an arbitrary time (normalized 
again to i = in this section) between the inception of the contract and its 
maturity T. As usual, it is obtained as difference between the default and fee 
legs, but, contrary to the usual pricing formulas, we take into account exactly 
the effect of the counterparty default leg. 

Notice that we insist in deriving exact formulas but, for practical purposes, 
it may be convenient to use the standard approximation schemes to simplify 
the use of the formulas and fasten the computations. For example, in order to 
compute the fee leg, it may be convenient to assume that defaults can occur 
only at times ti,i = l...n, so that the following integral expressions can be 
safely replaced by finite sums involving only the computation of the default 
probabilities P(ti A r2 <t). 

Let us write Ds for the "standard" default leg of the CDS, that is, the present 
value of the cash-flows corresponding to payments by the seller of protection if 
there is a default occurring on the assets underlying the CDS contract. The 
total default leg Dtot of the CDS contract is given by the sum Dg + Dc, where 
Dc is given by Thm. 



Lemma 3.4 The standard default leg is given by: 
Ds = C{l-Ru) / e 
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The proof follows from the same reasoning as the one of Lemma l3.2l and can 
be omitted. 
Notice that: 

Q(ri G dt, ri = n A ra, Zi{t) + iZ2{t) = d^i ■ e'") = 
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Let us write s for the spread of the CDS contract (the continuous coupon 
rate served by the protection buyer to the protection seUcr, tiU the end of the 
contract, or the default of the underlying, or the default of the protection seller); 
s is fixed at inception of the contract and remains constant till the maturity T. 
Recall that the fee leg of a credit derivative contract is the present value of 
cumulated payments by the protection buyer. We write from now on r for 

Lemma 3.5 The fee leg F of the CDS contract is given by: 

1 — p~^'^ 1 — p^^^ 

Fs-' = P(t > T) + E[ lr<T] 

r r ~ 
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The first term was computed in [T31 p. 698], the last two follow from our 
previous computations. 

Theorem 3.6 The fair price of a CDS, taking into account the counterparty 
risk, is given by: 

Dtot - F 

where F and Dtot Cife given by the previous Lemmas. 

4 First-to-default on two underlyings 

As was mentioned in the introduction, most credit derivative transactions on 
two default-prone instruments with time-dependent payoffs can be priced using 
the techniques developed in the present article. As an example, we derive the 
fair spread of a first-to-default contract on two underlyings. 

Most notations are as above, when we were dealing with the counterparty 
risk on a CDS, and we do not recall them excepted when necessary. The only 
difference is that, now, the two default-prone entities are treated on the same 
footing. We still write C for the notional value of the contract and assume a 
given recovery rate R (e.g. 40%) so that, at the first default occuring before the 
maturity of the contract, the buyer of protection will receive C(l — R) (recall 



this is the very definition of a first-to-default contract; on the other hand, the 
buyer of protection pays a continuous fee -still called the "spread"- till the first 
default occurs or till the maturity of the contract if no default occurs before the 
maturity) . 

Lemma 4.1 The default leg of the contract is given by: 
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On the other hand, the fee leg is given by the same formula as in Lemma l3.5l 
-this is because the coupon is served by the protection buyer till a first default 
occur in both cases (a CDS with counterparty risk or a first-to-dcfault swap on 
two under lyings). 

Theorem 4.2 The fair spread s of a fir st-to- default swap on two underlyings 
at inception of the contract is given by 

D 

where D is given by Lemma \4.1\ and F' = Fs^^ by Lemma \3.5[ 

5 Appendix 1: Theoretical default legs for CDS 

We include in this section a proof of the formulas for the market price of a CDS 
conditional to the default of the counterparty. These are classical computations 
and/or variations thereof. Wc give therefore only short indications and refer to 
[U Chap. 3] for further details. 

The notations arc as in Section 3; N stands as usual for the cumulated 
standard normal distribution. 

Lemma 5.1 On the set T2 < ti, we have: 



O"! \^T — Ti CTi -JT — T2 

where Yt := In(^) = In(^) + v^t + G^B^it), z^i := r - k^ - 71 - \al 
The lemma is a direct application of Cor. 3.1.1 in |4]. 
Lemma 5.2 For a,b,c G R with b < and c^ > a: 

r oT,,r/& — ex, d + c , , d~c,, . 

with d = ^/c? - 2a, g{y) := e''^''~'^^N{^), h{y) := e''('=+'^)iV(^). 



See Lemma 3.2.1, Fla. 3.16 in g]. 
Corollary 5.3 We have, on ti > T2: 

\/s-T2 \/s-T2 



where a := y^^ + 2r, f3 :== ^, /x^, := ^. 
Indeed by the previous lemmas. 

r2 CriVs-T2 CriVs-T2 

2a \/T - T2 2a ^T - T2 



_^g-2i'l<T 



2a \fT~-T2 2a ^/T — r2 



and the formula follows. 

6 Appendix 2: First hitting time in a polyhedral 
domain 

Our results rely largely on the following two-dimensional extension of a Theorem 
of Daniels, see [H Fla 3.1]. This natural extension is important in view of 
applications to multidimensional Brownian processes, particularly in the setting 
of credit risk. It was first announced without a proof in [TT] and could be easily 
extended to higher dimensional cases. We could not find a proof in the literature, 
and include therefore a short demonstration in this appendix. 

Let z^° (abbreviated to Zt when no confusion can arise), a planar Brown- 
ian starting at (xq, j/o) evolving in a polyhedral domain D with absorbing bound- 
ary dD. We write r for the first hitting time of the boundary and /(xq, j/O; 2;, y, t) 
(abbreviated to /(a;, y, t) when no confusion can arise) for the density of the sur- 
viving process: 

f{x,y,t)dxdy := P{Zt G {dx,dy),T > t) 

We look for the distribution of hitting times: P{t G dt,Zr G dp), dp G dD 
and can assume (up to a planar rotation) that dD is horizontal and an upper 
boundary for D in the neighborhood of p, so that dp is a horizontal line element: 
dp = {da, b), with p = (a, b) G dD. 

Theorem 6.1 We have: 



1 r\ 

P{T(^dt,Zr^{da,b)) = --—f{a,b,t)dadt. (1) 
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Sketch of the proof : Let us write (f>{x,y, At)dadt for P{t E dt,Zr E 
{da,b)\Zt-At = {x,y)), where in practice we will choose At << 1 and study 
the asymptotics when At goes to 0. From Daniels' theorem [SJ Fla 3.1] which 
computes the first exit density for a one-dimensional BM, together with the 
independence of the horizontal and vertical components of a two-dimensional 
BM, we get, in a neighborhood of (a, b): 

1 (a-x)2 (b-yf „ 

Then, 

P{t e dt, Zr e {da, b))da^^dt^^ = / /(.t, y, t - At)(j){x, y, At)dxdy 

D 

[/(a, y, t-At)+{x-a)f'^{a, y, t-At)+ ^^~"^ /^'(a, y, t-At)]c^{x, y, At)dxdy 

D 

By symmetry, in a neighborhood of (a, b) in D, (j){x, y, At) = (j){2a — x, y, At), 
so that the second term of the expansion vanishes. Since dD is an absorbing 
boundary, /(a, 6, t - At) = f^{a, b, t - At) = f^{a, b,t- At) ^0 and the third 
term reads: 

^ ^ f;^{a,y,t-At)cf,{x,y,At)dxdy ^ -— ^ ^ ' _--^ e-^^rfx 



2 ^"^ '"' '"^ '"' ' '^ Aty_^ 2 V2^Ai 

/■ft 



* / 6 - y 2-===e-^A^ ^ a, 6, < ~ At dy 
J-oo V27rAt ay 

0(A<) 



Similarly, the first term reads: 

/"ft iU \ 2 /"CXD -1 2 

/(a, y, f - At)</.(x, y, At)da;dy ~ / -^^El^e-^^^^ f{a,y,t - At)dy ^^^e-^^'^^dx 



D 



-)=-f^e-^^f{a,y,t~At)dy 

_oo y iTi At-^ I '^ 

C 4^e-'^fUa,b,t-At)dy 



2' 



= --f'yia,b,t~At) 



The Theorem follows. 
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